library(dplyr)
library(knitr)
library(maptools)
library(TreeSegmentation)
library(ggplot2)
library(rgl)
library(clue)
library(lidR)
knit_hooks$set(webgl = hook_webgl)
opts_chunk$set(warning=F,message=F)
#set color ramp for treeID
col = pastel.colors(200)
Load in ground-truth
shps<-list.files("/Users/ben/Dropbox/Weecology/ECODSEdataset/Task1/ITC/",pattern=".shp",full.names = T)
itcs<-lapply(shps,readShapePoly)
names(itcs)<-sapply(itcs,function(x){
id<-unique(x$Plot_ID)
})
itcs<-lapply(itcs,function(x){
proj4string(x)<-CRS("+init=epsg:32617")
return(x)
})
Example Pipeline
Read in Data
ground_truth<-itcs[[20]]
fname<-get_tile_filname(ground_truth)
tile<-readLAS(paste("/Users/ben/Dropbox/Weecology/NEON/cropped_",fname,sep=""))
tile@crs<-CRS("+init=epsg:32617")
plot(tile)